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ABSTRACT 

Only a few binary systems with compact objects display TeV emission. The physical properties of the 
companion stars represent basic input to understand the physical mechanisms behind the particle acceleration, 
emission, and absorption processes in these so-called gamma-ray binaries. Here we present high-resolution and 
high signal-to-noise optical spectra of LS 2883, the Be star forming a gamma-ray binary with the young non- 
accreting pulsar PSR B 1259-63, showing it to rotate faster and be significantly earlier and more luminous than 
previously thought. Analysis of the interstellar lines suggest that the system is located at the same distance 
as (and thus is likely a member of) Cen OBI. Taking the distance to the association, d = 2.3 kpc, and a 
color excess of E(B-V) = 0.85 for LS 2883, results in My ~ -4.4. Because of fast rotation, LS 2883 is 
oblate (Req ~ 9.1 Rq and Rpoie ~ 8.1 and presents a temperature gradient (Tiq « 27500 K, loggeq = 3.7; 
Tpoie ~ 34000 K, loggpoie = 4.1). If the star did not rotate, it would have parameters corresponding to a late 
O-type star We estimate its luminosity at log{L^,/LQ) ~ 4.79, and its mass at « 3OM0. The mass function 
then implies an inclination of the binary system /orb ~ 23°, slightly smaller than previous estimates. We discuss 
the implications of these new astrophysical parameters of LS 2883 for the production of high energy and 
very high energy gamma rays in the PSR B1259-63/LS 2883 gamma-ray binary system. In particular, the 
stellar properties are very important for prediction of the line-like bulk Comptonization component from the 
unshocked ultra-relativistic pulsar wind. 

Subject headings: binaries: close — gamma rays: stars — pulsars: individual (PSR B1259-63) — stars: 
emission-line. Be — stars: individual (LS 2883) — X-rays: binaries 



1. INTRODUCTION 

Three binary systems containing a massive star and a com- 
pact object that clearly display TeV emission are currently 
known: LS 5039, LS I H-61 303 and PSR B1259-63 (see, e.g., 
iRibo 2008). In PSR B1259-63 the compact object is a neu- 
tron star, first discovered as a 47.7 ms radio pulsar, in a very 
wide and eccentric o rbit (e = 0.87, Pmu = 1236. 9 d) around 
the Be stai- LS 28830 (iJohnston et al.lll99llT99l . A double- 
peaked non-thermal and unpulse d radio outburst take s place 
around each periastron passage (Johns ton et al.ll2005l) . This 
radio emission has recently been resolved at AU scales by 
iMoldon et al.l ( 1201 Ih . who place it outside the binary system. 
Double-peaked outbursts have also been detected in s oft and 
hard X rays, and at TeV energies by HESS ( Uchivam a et al.l 
120091: IGrove et al.l ll995tlAharonian et al.ll200 5. 2009). Upper 
limits were obtained at GeV energies by EG RET around the 
1994 periastron passage (iTavani et al.ll 19961) . 

In the currently preferred scenario, the cold ultra-relativistic 
wind of the pulsar interacts with the massive star wind within 



Partially based on observations collected at the European Southern Ob- 
servatory, Paranal, Chile (ESO 282.D-5081) and the South African Astro- 
nomical Observatory. 

- Star 2883 in the catalog of Luminous Stars in the Southern Milky Way 
(Stephenson & Sanduleak 1971) is also known as CPD -63°2495. The use of 
SS 2883 should be avoided, as this acronym refers to the catalog of emission- 
line stars by Stephenson & Sanduleak ( 1977), which contains only 455 stars. 



the binary system. Relativistic electrons are accelerated 
in a shock region where the pressures of both winds bal- 
ance. These electrons produce synchrotron radiation and 
upscatter UV-optical photons from the companion star via 
Inverse Compton (IC), giving rise to broadband emission 
up to TeV ga mma rays (e.g., Kirket al. 1999; Dubus 200^ 
iKhangulvan et alj I2007t lAharonian et all |2009). GeV and 
TeV emission can also be produced from IC scattering of stel- 
lar photons by electrons in the unshocked pulsar wind, de- 
pending on its Lorentz factor Radio emission outside the 
binary system is produced by long-lived particles traveling 
away in a kind of cometary tail. 

In these models, the companion star is the source of seed 
photons for IC scattering, expected to be the dominant ra- 
diation mechanism. Fundamental parameters, like the lumi- 
nosity or the effective temperature, have to be known ac- 
curately to properly model, and thus be able to understand, 
how this gamma-ray b inary p roduces the obser ved multi wave- 
length properties (see Khangulvan et al. 2007). LS 2883 was 
observed in the optical by Johnston e t al.. (1994), who identi- 
fied emission lines typical of a Be star. Based on the presence 
of strong He I emission, they concluded that it was earlier than 
B3. Ar guing that the star w as likely located in the Sagittarius 
Arm, Johns ton et al.l(ll994l) assumed a distance of 1.5 kpc and 
a spectral type B2e. 

Here we report on new astrophysical parameters of 
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LS 2883, a new distance estimate and a new inclination of 
the binary orbit. These resuhs will allow for better modeling 
of the multiwavelength data acquired close to the periastron 
passage of PSR B 1259-63 in 2010 mid December, the first 
covered by the Fermi and AGILE gamma-ray missions. 

2. OBSERVATIONS AND SPECTRUM DESCRIPTION 

Observations of LS 2883 were obtained on 2009 August 
28 usin g the Ultraviolet and Visual Echelle Spectrograph 
(UVES: lDekker et al. '2000) mounted on the 8.2 m Veiy Lai-ge 
Telescope (VLT) UT2 Kueyen at Cerro Paranal, Chile. UVES 
was used in dichroic 2 mode with cross-disperser CD#2 in 
the blue arm, giving coverage of the 373-499 nm range, and 
CD#4 in the red arm, allowing coverage of 565-946 nm with 
a small gap around 760 nm. We used Image Slicer #2 and a 
0"5 slit, giving a resolving power R ^ 80000. 

In addition, intermediate-resolution spectra of LS 2883 and 
three nearby blue stars were obtained on 2006 May 8, us- 
ing the unit spectrograph on the 1.9-m telescope at the South 
African Astronomical Observatory (SAAO) in SutherlandQ 
with grating #6 (600 lines mm"'). This configuration covers 
380-560 nm with R « 2000 (measured on arc frames). 

The spectrum of LS 2883 is t ypical of a Be star of e arly type 
with a well-developed disk (cf. iJohnston et al.lll994l) . Balmer 
lines have strong double-peaked emission components, with 
the red peak clearly stronger than the blue one (the two peaks 
are blended in Ha). Representative emission lines are dis- 
played in Figure [T] All Balmer lines higher than H7 clearly 
show the photospheric absorption wings, as the emission fea- 
ture is definitely narrower than the absorption line, allowing 
the determination of effective gravity. Paschen lines between 
Pa 11 and Pa 23 do not show any sign of the photospheric 
component. Other emission lines typical of early-B stars, such 
as O I 8448A or He I 5875, 6678, 7065A appear as double- 
peaked features. The 1 1115 A triplet is blended into a strong 
single feature. There are many weak emission lines, mostly 
from Fe II transitions, but also including Si II 4128 & 6347 A. 
The SAAO spectrum also shows strong emission lines from 
Fell 5019, 5169 &5316A. 

The equivalent width of Ha, EW(Ha), is -54 ± 2A, signifi- 
cantly higher than measured bv IJohnston et al.l (119941) . -40A, 
though the overall shape is very similar. The EW(H/3) is 
-4.5 ± 0.4 in the 2006 SAAO spectrum and -4.6 ± 0.2 in 
the 2009 VLT/UVES spectrum, showing no significant vari- 
ability. As seen in Figure [T] all the upper Paschen and 
Balmer features have similar peak separation, with values 
170±10kms-i. 

3. ASTROPHYSICAL PARAMETERS 

3.1. Distance and color excess estimates 

Since model fitting can only give values for T^ff and \ogg, 
we need an accurate distance to LS 2883 to derive M, and L,. 
We use interstellar atomic lines in the spectrum of LS 2883 
to study the radial velocity distribution of material along this 
Une of sight {£ = 304 .°2, b = -l.°0). Figure |2] displays three 
strong interstellar lines in the spectrum of LS 2883, the two 
components of the Na I doublet and one of the members of the 
Ca II doublet. All three lines (and other interstellar lines, like 
CN 3875A) show the same morphology, with two clearly sep- 
arated components. The first component, with low positive 

^ |http : //www ■ saao ■ ac ■ za/f acilities/instrumentation/ 
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Fig. 1. — Representative emission lines in the VLT/UVES spectrum of 
LS 2883, in velocity space. Velocities are measured with respect to the Local 
Standard of Rest (LSR) assuming Solar motion +16.6 km s"' towards Galac- 
tic coordinates £ = 53°, b = +25°. From top to bottom and vertically shifted 
for clarity, we show H»j (green), US (red), H/3 (black), and Pa 11 (blue). Their 
central wavelengths are quoted in A to the right of the figure. The dotted lines 
mark the position of the emission peaks in Pa 11, with velocities in kms"' 
annotated. 

velocity, is likely due to nearby clouds associated with the 
Southern Coalsack, which sho ws LSR veloc ities in the range 
-10 to -1-8 kms"' in CO maps CNyman et alJ ll989). The sec- 
ond feature is a combination of several weaker components 
with negative velocities (between -14kms"' and -36 kms"'), 
values produced by clouds in the Sagittarius Arm. Very 
similar components are observed in the interstellar lines of 
HD 112272 (£ = 303.°5, b = -l.°5) a B0.5Ia supergia nt be- 
lieved to be a member of Cen OB 1 (iHunter et al.ll2006l) . 

In this direction, the H I shell GSH 305-1-01-24, with 
dynamical d = 2.2 ± 0.6 kpc, has an average velocity 
-24 kms"', and is believed to be associated to Cen OBI 
jMcClure-Griffi^ths et aT.i 12002). LS 2883 shows interstellar 
components more negative than this value, suggesting that it 
is beyond this shell. In fact, the Galactic rotation curve would 
place LS 2883 at >2.5 kpc (see inset of Figure |2]i. More 
distant stars along this sightline show narrow components at 
« -50 kms"', believed to arise from clouds in the Norma- 
Centaurus Arm, at d >4 kpc ( Kaper et al.. 2006, and refer- 
ences therein). The lack of this component in LS 2883 sets an 
upper limit to its distance and firmly places it near Cen OBI. 

Photomet ric measurem ents of LS 2883 in th e ht- 
erature (Klare & Neckel 1977; Schild et alj IT981 
iWester lund & Garnier 1989; Drilling 1991) suggest some 
variability, at the level of a few hundredths of a magnitude, 
typical of Be stars. Taking (B-V) = 0.73 as representa- 
tive, the intrinsic color corresponding to the model fit (see 
below), (B-y)o =-0.28, implies E{B-V) = 1.01. Not all 
2£=^tt3ffig£wRI'f;ning is interstellar, as the disks of Be stars give 
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Fig. 2.— Interstellar lines in the VLT/UVES spectrum of LS 2883. The 
main panel shows the components of three of the strongest interstellar lines: 
Na I 5890A (solid black), Na I 5896A (dotted red) and Ca II 3934A (dashed 
blue) in velocity space. Peaks in the absorption and the exti'eme edges are 
annotated. The inset shows the Galactic rotation curve (black crosses) along 
this line of sight {£ = 304°2, b = -l.°0) according to Brand & Blitz 1 1993), 
using a circular rotation velocity at the position of the Sun (Jqc = 8.5 kpc) of 
220 km s"' . The green dashed line indicates the velocity of the most negative 
peak measured towards LS 2883. 

rise to circumstellar red dening due to stro nger emission at 
longer wavelengths (e.g.. lDachs et al.lll988[). The cor r elatio n 
between EW(Ha) and E^\B-V) from iDachs etal] ([T988b . 
valid solely for isolated Be stars, predicts E'^'^(B-V) = 0.11, 
higher than for most Be stars. Higher values are observed 
in Be/X-ray binaries, but values > 0.3 mag are generally 
associated with transient events (e.g., iReig et aLl l2007h . 
Therefore it seems sensible to assume LS 2883 to have 
interstellar Eifi-V) = 0.8-0.9 mag. 

We have also estimated the color excess using four inter- 
stellar diffuse bands (DIBs) with EW> 100 mA in the spec- 
trum o f LS 2883: A5780, A5 797 (iHerbiglll993l) . A6202 and 
A6614 (iRawlings et alj|2003l) . A4430, is difficuh to measure. 
Though there is some dispersion, all support E{B-V) between 
0.8-0.9 mag. We shall thus assume E{B-V) = 0.85. Members 
ofCenO Bl nearLS 2883 have ^(B-y) between 0.7 and 1.1 
(iHumphreys 1978). 

The extinction law in this direction d oes not differ signif- 
icantly from the average R = 3.1 law (IWinkleri 11997b . We 
observed two luminous stars close to LS 2883 from SAAO. 
For LS 2888 (7' away), we derive a spectral type B0.2IIL 
Using photometry from the literature (Klare & Neckel 19771 
[^hild et al. 1983) and standard calibrat ions fo r intrinsic col- 
ors (Fitzger ald. 1970 ) and magnitudes (iTumerlll^Oi) . we es- 
timate E\B -V) = 0.83 and d = 2.9 ± 0.5 kpc. For LS 2882, 
only 30" away from LS 2883, we derive a spectral type Bl IV. 
Unfortunately, there is no accurate photometry published for 
this star. Its 2MASS (J-Ks) color, however, indicates that it 
is reddened by an amount very similar to LS 2888 The spec- 
tral types and magnitudes are fully compatible with the three 
stars being at similar distances. 

Therefore direct measurements and indirect evidence are 
consistent with the idea that LS 2883 has interstellar E(B- 
V) = 0.85 ±0.05 and is a member of Cen OBI. The distance 
to this association has been studied by different authors, with 
most recent values converging towar ds d = 2.3 kpc, though 
values up to t/ = 2 5 kpc are reported (iKaltcheva & GeorgievI 
1 1 9941: iHumphrevs" 1 978', and references therein). We will as- 
sume that LS 2883 is located at the distance to the association. 
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Fig. 3.— The rectified VLT/UVES blue spectrum of LS 2883 (solid black) 
and its best model fit (dashed red). To illustrate the range of acceptable val- 
ues, in the top panel we also show models calculated with Tiff = 32 000 K 
and logg = 3.60 (blue dashes) and 4.10 (green dashes). In the lower panel, 
we show models with logg = 3.80 and T^{{ = 31000 K (blue dashes) and 
Teff = 34000 K (green dashes). Features used for the analysis are marked 
with dashed lines. 



TABLE 1 

Fitted and derived apparent stellar parameters forLS 2883 . 

Parameter Value 

Teff (K) 32000^^^" 

logg 3.80+03« 

R,i (Rq) 9.0!| f 

log(L./L0)' 4.88!<]]« 

M.1(M0) 21.3!22j° 

Af«'.2(Me) 26.6!^^;^ 
"Assuming d = 2.3± 0.4 kpc and EiB-V) = 0.85 ± 0.05. 
*'Deproje cting the rotation al velocity with i = 35° , after applying the correc- 
tion of lFremat et all f200% . 

even though values as high as 2.8 kpc seem compatible with 
the data available. Values lower than 2.0 kpc seem to be ex- 
cluded by the kinematic data contained in the interstellar lines 
to LS 2883. 

3.2. Spectrum modeling 

Stellar parame ters have been calculated using Fastwind 
dPuls et al.l I2005h . a spherical non-LTE model atmosphere 
code with mass loss, developed to model detailed spectral line 
profiles of hot and luminous stars. The best fit model is shown 
in Figure [3] 

We determined the stellar rotational velocity using the 
Fourier transform method. We first degraded the spectrum 
to R= 10000 for a better SNR. Because of the lai-ge broaden- 
ing, many lines are blended. We took unblended line wings 
and mirrored them (thus producing a symmetric line) be- 
fore calculating the Fourier transform. This procedure gives 
results consistent with measurements for the blended lines. 
Our values, which rest on the metallic lines O II 4254A and 
C II 4267 A, are Vrotsin/ = 235 ± 15kms"' (projected rotational 
velocity; to be refined further below) and = 1 15 ± 50kms"' 
(so-called macro turbulence). Errors simply represent the dis- 
persion of individual values. 

Stellar parameters were obtained by visual fitting to 
the Balmer, He I/He II and Si Ill/Si IV Unes (see, e.g., 
iRepolust et allllOOllSimon-DfazllIoTol) . The fit is not as ac- 
curate as in these works, since the heavy veiling by emis- 
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sion lines makes the continuum definition a difficult task 
and some diagnostic lines either have emission components 
(such as He 1 4471 A) or are affected by nearby emission lines 
(as, for example. Si III 4552A). The best fit is obtained for 
Teff = 32000!f;jjgK and \ogg = 3.8^^1, where the uncertain- 
ties have been obtained following the procedures discussed 
in Sect. 6 of Repolust et al. (2004). A normal He abundance 
(A^(He)/(A^(H)-l-A^(He))= 0.09) is consistent with the observed 
spectrum, although a slightly increased He abundance is pos- 
sible. The lines typically used to estimate the mass-loss rate 
are strongly contaminated by emission components, prevent- 
ing us from obtaining an accurate value. We have adopted the 
mass-loss rate expected for a spherical star of the luminosity 
derived below according to the recipes of Vinketal. (200d|), 
resulting in a thin wind with no effect on the emergent pro- 
files. 

Using t/ = 2.3 ±0.4 kpc, £(B-y) = 0.85 ±0.05 and 7? = 
3 . 1 ± 0. 1 we derive My = -4.4 ± 0.4. From the model fit spec- 
trum and the absolute magnitude, we derive the values for 
R^,, L, and listed in Table[T] The fitted parameters corre- 
sp ond to an 09.5 V sta r according to the observational scale 
of iMartins et al.l (l2005h . in good agreement with the morpho- 
logical spectral classification that we derive from our spectra. 

However, for fast rotators, there are strong differences in 
effective gravity between the polar and equatorial regions, 
resulting in important temperature gradients and oblateness 
fTownsend et al.li2004 a nd therein). The bin ary mass func- 
tion /(Mns) = 1.53M0 jJohnston et alj|1994) implies an or- 
bital inclination angle of i'nrh = 24.°7;^ ^ '4 for a standard neutron 
star mass of 1.4M0. iMelatos et all (119951) find a tih of 10° 
between the orbital plane and circumstellar disk, expected to 
lie on the stellar equatorial plane. Therefore we are seeing 
LS 2883 under a low (^ 35°) inclination and the parameters 
obtained must mainly reflect the characteristics of the hotter 
polar region. 

For a centrally condensed star, the critical rotational angular 
velocity (gravity and centrifugal forces are equal at the equa- 
tor) is given by: 

ne. = ,/^. (1) 

A star rotating at f lmt will have R^q =1.5 Rpok 
(ICranmer & Owockil[T995h . In the case of LS 2883, with a 
deprojected Vrot = 410kms"', a rough estimation suggests that 
7?eq > l.l^poie and Lu = 51,/r2crit ^ 0.85, meaning that 2D ef- 
fects will be important. 

For a start, in fast rotators, both the Full Width at Half Max- 
imum (FWH M) of lines and the F ourier method underesti- 
mate Viotsin; (iTownsend et al.ll2004i Fremat et al. 2005). We 
therefore correct Viotsin/ using Fig. 8 of iFremat et al.i (i2005l) . 
finding v^^^sini = 260 ± 15kms"'. The deprojected velocity is 
then v^Q, « 450 km s"', even closer to the critical value. 

Though a complete 2D treatment of the stellar surface is 
beyond our scope, we can perform a quick evaluation of the 
effects of fast rotation. From the Teff and \ogg derived from 
our analysis (apparent values), we estimate the Teff, log^ and 
My of a non-rotating star that would give the same apparent 
parameters when rotating at the same Vi-ot as LS 2883. For 
th is, we use an itera tive method base d on the calculation s 
of IFremat etaTI (12005 ) and the tables of lCoUins etall (11991b . 
With these new values, we calculate other stellar parameters 



TABLE 2 

Estimated artwa/ stellar parameters forLS 2883. 



non-rotating 


rotatin 
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^eff 


33 500 K 


(' 


33° 




4.0 


UJ 


0.88 




-4.47 


Tiff (eq) 


27 500 K 


Rl 


9.2R0 


Tiff (pole) 


34 000 K 


log(LO/LQ) 


4.98 


\ogg (eq) 


3.7 


M» 


31 Mq 


\ogg (pole) 


4.1 






Req 


9.7 Rq 






Rpole 


8.1 Rq 






\og(U ILq) 


4.79 



and iterate again to take into account the change in the inclina- 
tion due to the change in . The results of this procedure are 
the stellar parameters that LS 2883 would have if it did not ro- 
tate (Table|2] left panel) and its most likely actual parameters, 
taking into account fast rotation (Table |2] right panel). 

LS 2883 presents important differences in T^a and log^ be- 
tween the polar and equatorial regions, rendering the apparent 
parameters little more than an approximate guess. The hypo- 
thetical non-rotating LS 2883 would have parameters roughly 
con^esponding to an 08V star. Its observed spectrum is a 
consequence of fast rotation. The mass derived, 31M0, is 
somewhat high for these parameters, but is subject to large 
uncertainties, as it depends strongly on the corrections for fast 
rotation and the distance assumed. 

4. DISCUSSION 

Our spectroscopic observations and atmosphere model fit- 
ting have provided new physical parameters for LS 2883, the 
massive star forming a gamma-ray binary with the young non- 
accreting pulsar PSR B 1259-63. The higher temperature and 
luminosity of the optical star presented here, Teq ~ 27 500 K, 
Tpoie ~ 34000 K, U = 2.3 x lO^** erg s"', as compared to 
previous estimat es, T^s = 23 000-2 7 000 K, U = (0.3-2.2) x 
10^^ erg s"' (see lKhangulvan et alj[2007 ). imply a significant 
revision of the parameters and conditions for the production of 
nonthermal radiation in this binary system, especially in the 
gamma-ray band. High (GeV) and very high (TeV) energy 
gamma-rays from this system can be produced in two distinct 
regions: a) in the unshocked pulsar wind, i.e., a cold ultra- 
relativistic outflow expanding outwards from the pulsar with 
Lorentz factor F < 10^; and (b) in the region of the terminated 
pulsar wind. In both regions, the dominant gamma-radiation 
mechanism is IC scattering of relativistic electrons. While 
the Very High Energy (VHE) emission detected by HESS is 
most likely linked to the multi-TeV electrons accelerated af- 
ter termination of the wind, GeV gamma-ray emission can be 
effectively produced also by the unshocked pulsar wind. The 
higher luminosity of the optical star obviously implies an en- 
hanced interaction rate and consequently higher luminosity of 
IC gamma-rays {L^ oc L*). The Fermi AGILE gamma-ray 
missions should be able to detect, at epochs close to the pe- 
riastron, the line-type emission of the unshocked pulsa^r win d 
(iKirk et alJll999l:lBall & Ki^lIOOOHKhangulvan et alJ[2007h . 
and in this way measure its Lorentz factor and mechanical 
power 

At first glance, the enhanced luminosity of the optical star 
should have a weaker impact on the gamma-ray emission re- 
lated to the termination of the pulsar wind, since gamma-ray 
production is expected to proceed in the saturation regime. 
However, the enhanced optical luminosity does affect the 
gamma-ray emission for two reasons. Firstly, it introduces 
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a significant orbital dependence of the gamma-ray flux due to 
the Compton deceleration of the unshocked pulsar wind. This 
effect has been realized bv'Khangulvan et al] ( 1200 7). who in- 
dicated that the deficit of VHE emission observed by HESS 
close to the periastron passage may be explained by this ef- 
fect, provided that the luminosity of the optical star is as large 
as '--^ 4 X 10^** erg s~' . Although the value found here is smaller 
by a factor 1 .7, the Compton drag effect should be combined 
with gamma-gamma absorption, since for the updated orbital 
inclination the maximal attenuation (by a factor of 2) occurs 
close to periastron passage. We note, however, that given 
the large orbital separation the total absorbed energy will be 
small, thus the contribution from the electromagnetic cascade 
should be negligible. 

Finally, the production rate of VHE gamma-rays is affected 
by the increase of the stellar temperature and the change of the 
orbital inclination. The higher photon temperature leads to a 
more pronounced impact of the Klein-Nishina effect. This, 
together with the smaller orbital inclination obtained here, 
should result in a weaker orbital phase dependence of the 



VHE gamma -ray production, in agreement with recent HESS 
observations dAharonian et al.ll2009l:lKersclihagglll201 1 ). 

It is clear that the new physical parameters of LS 2883 re- 
ported here have to be taken into account in the interpretation 
of the high and very high energy gamma-ray observations of 
the gamma-ray binary LS 2883/PSR B 1259-63, such as the 
multiwavelength campaigns close to the 2010 December pe- 
riastron passage, which include for the first time the partici- 
pation of the Fermi and AGILE gamma-ray missions. 
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